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ABSTRACT 

A substantial population of high redshift early-type galaxies is detected in very deep UB- 
VRIJHK images towards the HDF-South. Four elliptical profile galaxies are identified in the 
redshift range l<z<2, with very red SEDs, implying ages of ^2 Gyrs for standard passive evo- 
lution. We also find later type IR-luminous galaxies at similarly high redshift, (10 objects with 
z>l, H<25), with weak UV emission implying single burst ages of ^1 Gyr. The number and 
luminosity-densities of these galaxies are comparable with the local E/SO-Sbc populations for 
Q„^>0.2, suggesting that the major fraction of luminous Hubble-sequence galaxies have evolved 
little since z ~ 2. A highly complete photometric redshift distribution is constructed to H=25 
(69 galaxies) showing a broad spread of redshift, peaking at z~1.5, in reasonable agreement 
with some analyses of the HDF. Four 'dropout' galaxies are detected at z ~ 3.8, which are 
compact in the IR, ~0.5 kpc/h at rest 3500A. No example of a blue IR luminous elliptical is 
found, restricting the star-formation epoch of ellipticals to z>10 for a standard IMF and modest 
extinction. 

Subject headings: galaxies: evolution 



1. INTRODUCTION 

The detection of early type galaxies at z>l requires 
deep IR imaging to overcome the observed large rest- 
frame optical K-correction and correspondingly mini- 
mal inferred passive spectral evolution for this class of 
galaxy (e.g. Stanford, Eisenhardt & Dickinson 1998). 
The highest redshift example of an early-type galaxy 
is a luminous object reported by Dunlop et al. (1996) 
at z=1.55. Its restframe UV spectrum is dominated 
by late-type F-stars, equivalent to an age of « 3Gyrs, 
for standard stellar synthesis models (Spinrad et al. 
1997). Following ordinary ellipticals to z > 1 requires 
much deeper IR imaging now possible from space. 
The NICMOS observation in the Hubble Deep Field 
South (Fruchter et al 1999) represents a three magni- 
tude increase in depth over the deepest ground based 
IR data, and with usefully high spatial resolution. 
This, in combination with deep optical images from 
the VLT allows easy identification of red L* ellipti- 
cals to z ~ 2.5, where significant spectral and density 
evolution is widely anticipated. 

In section §1 we briefly outline our spectral and pro- 



file fitting procedures and simulations as applied to 
the NICMOS and VLT observations. In §2 we discuss 
detections of early-type and other IR bright galaxies 
at high-redshift. In §3 we compare their statistical 
properties with local populations of Hubble sequence 
galaxies of similar spectral type, for different evolu- 
tionary models. In §4 we briefly discuss a the absence 
of bright blue precursor ellipticals in the IR. 

2. PHOTOMETRY, REDSHIFT ESTIMATION AND 
PROFILE FITTING 

Three sets of images are utilized, the HDF-S NIC- 
MOS observations in JHK (Fruchter et al. 1999), the 
higher resolution unfiltered HST STIS image (Gard- 
ner et al. 1999) and the VLT Test Camera images in 
UBVRI (Fontana et al. 1999). We use the NICMOS 
pipeline zero-points and have computed the VLT ef- 
fective zero-points for the combined optical images 
from the ESO photometric solution. The filter func- 
tions and instrument response curves are kindly made 
available by ESO and STScI. The detection is per- 
formed using SExtractor (Bertin & Arnauts 1996) 
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in a combined J + H frame. Near total magni- 
tudes are measured in J,H,K and PSF matched STIS 
images within a matched aperture set. A total of 

200 galaxies are detected to H=27 within the cen- 
tral 50" x 50" field. The VLT UBVRI images have 
PSF's ~ 3 times wider that the NICMOS exposures, 
and require a different procedure. We measure mag- 
nitudes within a / /9 aperture and correct them to 
near-total magnitudes using versions of the J + H 
band image degraded to the PSF of each of the VLT 
passbands to ensure the colors are defined within the 
same apertures for all bands. 

Photometric redshifts are estimated using the 
methods described in Bemtez (1999). The empirical 
template set of Coleman, Wu, & Weedman (1980) 
is used, and augmented with bluer local starbursts 
(Kinney et al. 1996), to determine the redshift prob- 
ability distribution for each galaxy with Bayesian 
weighting. The fraction of objects with highly reli- 
able photometric redshifts (p>0.99) is 57/69 to H<25 
(24/25 for objects classified as earlier than Scd). 
Fainter than this magnitude limit the reliability of 
the photometric redshift estimates quickly degrades. 
We also allow for the quoted 10% uncertainty in the 
NICMOS zero-points for our spectral fits. 

For redder early types, more detailed fits to single 
burst synthetic models are justified to accommodate 
the known spread of metallicity and age. For this, the 
compilation of Leitherer et al. (1996) is used, includ- 
ing the synthetic galaxy spectra of Bruzual & Chariot 
(1993). We also check that the colors of these objects 
are not compatible with dust reddened star-forming 
galaxies by generating reddened versions of the late 
and Im types in our template set. 

Profile fitting is performed for bright galaxies in 
the IR (H<23.7). For this we simulate NICMOS im- 
ages using the methods described in Bouwens, Broad- 
hurst &: Silk (1998) to account for the pixelization, 
dithering and point spread functions and the avail- 
able weight map for the NICMOS dithering pattern. 

3. HIGH- REDSHIFT RED GALAXIES 

The H < 25 sample contains five clear elliptical 
galaxies, of which four lie above z=l. Their SED's 
and surface brightness profiles are compared with the 
models in Fig 1 and listed in Table 1. All the ob- 
jects spectrally classified as ellipticals are also found 
to have spheroidal profiles (Fig 1). The most impres- 
sive object (NIC3/ET4, Fig 1, Table 1) is extremely 
red, consistent with minimal passive evolution at an 
estimated z = 1.94±0.15, with a clear de-Vaucouleurs 
profile and uniform color over its surface (~ 300 pix- 
els). The effective radius 1.6(2.2)kpc/h is consistent 
with its inferred luminosity Mab = — 21.96(— 22.6) 
in restframe R for Q m = 1(0.1) (Binggelli, Sandage 



& Tarenghi 1984). This IR conspicuous object is ab- 
sent or only marginally detected in the U,B,V,R,I and 
STIS bands (Fig 1). Dust reddened spectra extin- 
guished with standard reddening (Calzetti, Kinney, 
& Storchi-Bergmann 1994) are unable to generate the 
steep break from the relative flat IR. Another ellip- 
tical identified at z = 1.66 ± 0.25 (NIC3/ET3, Fig 
1, Table 1) with a clear de-Vaucouleurs profile, is the 
luminous object reported by Treu et al. (1998) in 
the earlier NICMOS test images and confirmed by 
Stiavelli et al. (1999). 

We also find several IR-bright high-redshift galax- 
ies of later spectral type, 10 objects with z>l and 
H<25. Two typical examples are shown in Fig 1. 
The profiles of these objects are generally exponen- 
tial. The SED's are very conspicuous with a turn up 
towards the UV, and a distinctive bump in the IR 
separated by a trough, so there is little uncertainty 
in redshift (e.g. Fig 1). The best fit single burst 
minimum ages are typically ~ 1 Gyr. Comparison 
with the SEDs of CWW show these objects span the 
spectral types of local Sab to Sbc galaxies. The full 
sample of 24 galaxies with H<25 (open histogram, 
Fig 2) have a redshift distribution peaked at z 1.5 
with high-z tail similar to that found for the HDF 
(Connolly et al. 1997, Fernandez-Soto, Lanzetta & 
Yahil 1998, Bemtez 1999). 

4. DENSITY EVOLUTION 

A reliable measurement of the space density of mas- 
sive galaxies at high redshift provides a basic check 
of hierarchical models of galaxy formation, where the 
most massive haloes form late (Lacey Sz Cole 1993) 
and in which ellipticals and massive spheroids may be 
viewed as the end product of the merging of smaller 
systems. 

The presence of so many luminous red galaxies at 
high-z in such a small field is surprising and indi- 
cates that relatively evolved luminous galaxies are 
commonplace at z>l. Using the Pozzetti, Bruzual, & 
Zamorani (1996) local luminosity functions for E/SO- 
Sbc as input, we recover 4(9) ellipticals from our re- 
alistic NICMOS simulations for z > 1, H < 25, with 
Q m = 1(0.2) and no evolution. Adding the later- 
type IR-luminous galaxies allows a useful comparison 
with models (Fig 2). The redshifts and luminosities 
observed above z=l are consistent with no density 
or luminosity evolution of standard Hubble sequence 
galaxies for 0.2 < VL m < 1. Note, passive evolution is 
minimized here for consistency with the red SEDs by 
choosing Hq = 50km/s/Mpc and Zf = 15. The hier- 
archical model parametrization of Kauffmann, Char- 
lot & White (1996) for f2 = 1 underpredicts the num- 
bers as expected. Hierarchical models with Q, < 1 
may allow a declining density but must be checked 
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for self consistency. The optical- IR colors prefer older 
higher redshift models and are consistent with no size 
evolution, but do not distinguish clearly between the 
above cosmological models. 

The integrated luminosity density to H=25 at 
l<z<2.5 is comparable with local estimates inte- 
grated over E/SO-Sbc types (see Fig 2). The 
predicted luminosity density for z > 1 at rest- 
frame 6000A is 4.3 x 10 19 ergs/cm 2 /s/Hz compared 
to 3.9 ± 1.4 x 10 19 ergs/cm 2 /s/Hz, 2.41 ± 0.9 x 10 19 
ergs/cm 2 /s/Hz and 2.06 ±0.8 x 10 19 ergs/cm 2 /s/Hz, 
for n m = 1, n m = 0.1 and VL m + ft A = 0.25 + 0.75 
respectively. A factor of two decrease in luminos- 
ity density is implied relative to the popular low £l m 
models assuming passive or no stellar evolution. 

Our results seem to extend to z ~ 2 the results 
of other authors notably Lilly et al. (1995) and Ko- 
dama, Bower & Bell (1999) who find no evidence for a 
dramatic evolution in the density of early types galax- 
ies to z ~ 1 . The contrast of our results with the ab- 
sence of distant ellipticals claimed by others (Kauff- 
mann & Chariot 1998; Zepf 1997; Franceschini et al. 
1996; Barger et al. 1998; Menanteau et al. 1998) can 
be put in context by appreciating that the NICMOS 
data is complete to H « 27, which is 3 magnitudes 
fainter than the deepest ground based data and 6 
magnitudes fainter than most work on this question. 
Two of the four ellipticals found here in the IR would 
barely register in the I band (I~27) and are certainly 
not bright enough for morphological classification in 
the optical even from space. Recent claims of a de- 
cline in the elliptical density for z>l rely on signifi- 
cant passive evolution, where brightening enhances 
the expected numbers. The claim by Kauffmann, 
Chariot, k White (1996), based on the CFRS (Lilly 
et al. 1995) is subject to redshift incompleteness, in 
particular of red high-z galaxies, which are very hard 
to measure spectroscopically. Some of the red galax- 
ies reported here are bright enough in the IR to be 
detected in previous deep Keck images (Moustakas 
et al. 1997; Hogg et al. 1997) but morphological 
classification is not possible in ground based images 
and deep complementary optical data are required to 
establish the break. The NICMOS images directly 
show that luminous early type galaxies are present in 
significant numbers at high redshift with SEDs too 
red to be recognized as ellipticals in the optical even 
in the deepest WFPC2 images. 

An obvious problem with this kind of work is the 
presence of clustering in small field, rendering the re- 
sults uncertain. To explore this we have searched for 
spikes in the redshift distribution using the approach 
described in Bemtez (1999). The minimal number of 
structures along the light-of-sight to H=25 is three, 
at z = 0.55 (consistent with the 0.58 redshift spike 



found in the HDF-S by Glazebrook et al. 1999), 
z = 1.35 and z = 1.95. Each of these structures 
corresponds approximately to the redshift of one or 
more elliptical galaxies, and the two higher redshift 
peaks contain ss 1/3 of all the z > 1,H < 25 galax- 
ies. Our statistical analysis shows that it is virtu- 
ally impossible that all the early type galaxies at 
1 < z <; 2 belong to a single spike. A rough esti- 
mate of the expected clustering noise can be obtained 
by integrating the faint galaxy correlation function 
over the NICMOS field for the range l<z<2. Con- 
nolly, Szalay &; Brunner 1998 measure an amplitude 
A(10") = 0.12(a; oc #~ - 8 ) for the angular correla- 
tion of HDF-N galaxies within Az = 0.4 slices to 
z ~ 1.5. Assuming that A oc 1/Az, the total vari- 
ance in the number of galaxies in a redshift slice Az 
isa 2 (Az) ~ N(Az)[l + N(Az)(0A/Az)0.07\. For the 
interval 1 < z < 2 with 35 galaxies this yields a vari- 
ance of a = 8.3 galaxies, which is 40% greater than 
Poisson. If the early types are clustered like the inte- 
gral of all types then we estimate their number to be 
14 ± 4 in the range 1 < z < 2. 

5. VERY HIGH- REDSHIFT BLUE GALAXIES 

The lack of spectral evolution of early type galax- 
ies demonstrates that star-formation does not extend 
over time, but is confined to an early brief period, so 
that we should expect the precursors of these galaxies 
to appear extremely luminous and blue at higher red- 
shift. The width of this period is constrained to be 
short by the tight E/SO color-magnitude and M/L 
sequences (e.g. Bower, Lucey &; Ellis 1992, Van- 
Dokkum et al. 1998). More directly, the lack of blue 
ellipticals in the HDF (e.g. Franceschini et al. 1998) 
means this epoch lies at z > 5, and therefore may ap- 
pear in the near-IR for redshifts in range 5<z<15 as I 
or J band 'dropouts', provided the IMF of early-type 
galaxies is not unusually deficient in hot stars. 

Our simulated NICMOS images show that signif- 
icant numbers of elliptical galaxies are expected in 
the IR if they form luminous stars high redshift and 
were not significantly extinguished during the first 
Gyr. The simulations include an exponential burst at 
zj = 15 with an e-folding time of lGyr and predicts 
6, 10 and 15 objects for z>5 to H=26 for, Q m = 1, 
$7 = 0.1, Q m + Q\ = 0.25 + 0.75. In comparison, no 
I or J band 'dropout' galaxies are observed brighter 
than this limit. 

At lower redshift, five compact optically blue 
dropout galaxies are found, all consistent with a sin- 
gle redshift of z « 3.8. Fits to their spectra are shown 
in Fig 3. These objects have low SFR of 6-23M Q /yr 
and very compact sizes in the H-band, ~0.5 kpc/h 
at restframe 3500^4 (Q m =l), and therefore appear 
unlikely to be the precursors of luminous elliptical 
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galaxies if no evidence of SF is to remain by z~l. 

6. CONCLUSIONS 

We have detected a substantial population of early 
type galaxies at z>l. These include r 1 / 4 profile galax- 
ies in the redshift range z>l, all of which have sur- 
prisingly red colors despite their large redshifts. Later 
spectral-types are also present in this redshift range, 
with a conspicuous IR bump. These galaxies are 
all luminous and well evolved, with minimum single- 
burst ages >lGyr for a standard IMF. 

The clear detection of early type galaxies and 
spheroidal components of other high-z disk galaxies 
points to the existence of a roughly similar density 
of luminous Hubble sequence galaxies at z ~ 2, com- 
pared to the local Universe. The declining density of 
luminous galaxies expected in hierarchical schemes 
can be made consistent with the data only with for 



low £l m models. Low £l m helps accommodate the ab- 
sence of spectral evolution, but in any cosmology the 
lack of bright IR-blue precursors of elliptical galaxies 
requires a very high redshift of formation z>10, for a 
standard IMF and extinction. 

Although the field studied is small, the number of 
red high-z galaxies detected is surprising and strongly 
suggests that the bulk of luminous galaxies were 
formed at z ;> 2. In the short-term a deep NICMOS 
pointing in the HDF will provide similar data and 
help average over clustering. Progress on this ques- 
tion awaits deep diffraction-limited IR imaging on 
large telescopes and improved IR facilities in space. 

We would like to thank Bob Williams and the Hub- 
ble Deep Field South Team as well as Alvio Renzini 
and the VLT-UT1 Science Verification Team for mak- 
ing available this useful data. 
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Table 1 

z > 1 Early type galaxies in the HDF/NIC3 field 



name 


X 


Y 


mFimw,AB 


Zphot 




STISab - FWOWab 


Age(Gyr) 


NIC3/ET1 


181.6 


868.87 


21.24 ±0.01 


1.41 ±0.20 


0"34 


3.3 


> 1.5 


NIC3/ET2 


439.2 


352.34 


24.50 ± 0.06 


1.55 ±0.20 


0"31 


> 4.0 


> 1.5 


NIC3/ET3 


607.4 


371.52 


21.79 ±0.01 


1.66 ±0.25 


0"26 


4.2 


> 2 


NIC3/ET4 


646.8 


898.32 


23.18 ±0.02 


1.94 ±0.15 


0739 


> 5.4 


> 3.5 
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Fig. 1. — The left panel shows our best fit elliptical galaxies at z>l. On the left is the SED fit to the multiband photometry 
where the data points should be compared with the integrated flux over the passband from the best fit model and represented 
as open squares. The corresponding model spectrum is also shown. On the right are the convolved exponential (dotted) and 
de-Vaucouleurs (solid) best fits to the data, plotted as log surface-brightness versus linear radius for clarity. The two lower panels 
are representative examples of later spectral types at high redshift (Sab and Sbc), demonstrating the utility of U through K 
broadband information for identifying redshifts of these galaxy types at high redshifts. The positions of objects ET1 through ET4 
are listed on Tablel; objects SP1 and SP2 have coordinates (595.8, 1001.3) and (835.4, 1022.2). 
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Fig. 2. — Models of local E/SO-Sbc luminosity functions recovered from simulated images matched to the properties of the 
NICMOS South data. These are compared with the redshift distribution of E/SO-Sbc spectral types to H = 25 in the lower 
histogram and against apparent magnitude to the left. The predicted range of redshifts is similar to the data (shaded histogram) 
for all models, but the observed numbers are greater than expectations of hierarchical models for massive galaxies with Q m = 1 
and underpredict no-evolution or mild passive evolution with larger volume cosmologies. Consistency is found with mild or no 
stellar evolution for relatively flat Q m dominated models. Several Monte-Carlo realizations of the redshift-magnitude distribution 
for a no-evolution Q — 1 model (small points) are compared to the observations (big points). Dotted histograms give the redshift 
distribution for all types to H — 25. The lower panel shows the luminosity density with redshift for the spectrally early-type 
galaxies (solid squares) and for the whole sample (open squares) to H = 25 at rest 6000A Compared with the local luminosity- 
density of E^SO-Sbc galaxies (left most point), there is no obvious trend. The curve is our calculation of the luminosity density 
at rest 6000A using the Madau, Pozetti & Dickinson (1998) SF history for the Salpeter IMF and no dust. 
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Fig. 3. — The optical dropout galaxies showing the best single burst minimum ages ranging from 10 -10 yrs and star formation 
rates of between 6-23M©/yr for a Salpeter IMF and £1 = 1. These blue objects are very compact in the IR ~0.5 kpc/h at rest 
3500A and all consistent with lying at z ~ 3.8. No example of a higher redshift IR-dropout galaxy up to H < 26 is found, placing 
a useful constraint on single burst models of luminous galaxy formation (see text). D1,D2,D3 and D4 have respectively positions 
(244.3, 598.1), (318.3, 855.7), (857.6, 847.4) and (497.8, 770.7) 



